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Abstract 

Heating induced by an oscillating modulation of the interaction strength in an atomic Fermion 
pair condensate is analyzed. The coupled fermion-boson model, generalized by incorporating a 
time-dependent intermode coupling through a magnetic Feshbach resonance, is applied. The dy¬ 
namics is analytically characterized in a perturbative scheme. The results account for experimental 
findings which have uncovered a damped and delayed response of the condensate to the modula¬ 
tion. The delay is due to the variation of the quasiparticle energies and the subsequent relaxation 
of the condensate. The detected damping results from the excitations induced by a nonadiabatic 
modulation: for driving frequencies larger than twice the pairing gap, quasiparticles are generated, 
and, consequently, heating sets in. 

PACS numbers: 03.75Ss, 05.30.Fk 
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I. INTRODUCTION 


The realization of the crossover from a molecular Bose-Einstein condensate (BEC) to a 
Bardeen-Cooper-Schrieffer (BCS) superfluid of atom pairs jl-12| has opened the way to a 
variety of experiments on fundamental effects in quantum statistics and many-body physics. 
Essential to the versatility of this scenario has been the possibility of varying the interaction 
strength via a Feshbach resonance (FR), which has allowed the characterization of the effects 
under controlled conditions. In parallel, the need of explaining emergent phenomenology has 
brought about active theoretical work in the field. Recent research has dealt with the role 


of thermal fluctuations and the characterization of nonequilibrium situations 
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In this 


line, here, we extend previous theoretical work on the effects of an oscillating modulation 
of the interaction strength in a two-component Fermi gas of atoms {3- Our objective is 
to complete the analysis of the experiments of Ref. [15]. In them, a gas of ultracold ^Li 
atoms was prepared in the BCS regime through a magnetic FR, specifically, the (broad) 
FR at 834 G between the two lowest hyperfine states. The application of a sinusoidal 
modulation of the magnetic field was shown to lead the condensate fraction to oscillate with 
the driving frequency. Moreover, the oscillations were found to be damped and delayed 
with respect to the modulation, the damping time being much longer than the driving 


period. In previous work [3; the focus was put on the mechanism responsible for the 
delay: the modulation was shown to drive the system to an out-of-equilibrium situation, 
the deferred response being rooted in the finite relaxation time of the condensate. Here, we 
will concentrate on understanding the decay of the condensate fraction. To this end, we set 
up a framework where, through a partial analytical characterization of the dynamics, the 
origin of the damping processes can be identified. The decay of the oscillations will be linked 
to heating induced by the generation of quasiparticles. Nonadiabaticity, (on the gap time 
scale), will appear as a crucial component of the excitation process. Our scheme will allow 
us to trace the differential aspects of the mechanisms responsible for heating and delay in 
the system response. 


The outline o 
boson model 16 


re paper is as follows. In Section II, we generalize the coupled fermion- 


2l| by incorporating a time-dependent intermode detuning, or, equivalently. 


a time-dependent coupling. To deal with this variation of the basic model, a perturbative 
scheme based on the Hartree-Fock-Bogoliubov (HFB) description is developed. In Section 
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Ill, the relaxation of the condensate and the associated delayed response to the modulation 
are tackled. The heating effects are evaluated in Section IV. In Section V, our conclusions 
are summarized. 


II. THE FERMION-BOSON MODEL WITH A MODULATED COUPLING 
STRENGTH 


Our system consists of a gas of ultracold fermionic atoms with two hyperfine states 
coupled to a molecular two-particle state through a magnetic FR. To describe it, we apply 
the coupled fermion-boson model 


16 


18 


20|. The grand-canonical Hamiltonian reads 


H — ^ + hint ^ ®a+k,t®a-k,4.®f-k',t®f+k',t 


k,cr 


q,k,k' 


^ « + hvo) + gY^ (&q«Lk,t«a_k,i + h- 


( 1 ) 


q,k 


where /i is the chemical potential, N is the total number of atoms, (®k,cr) denotes 

a fermionic creation (annihilation) operator of an atom with momentum k and spin a, 
(a G andfol (6q) is a bosonic operator that creates (destroys) a molecule with momen¬ 

tum q. It is assumed that the two hyperfine states are equally populated. £k = h,^k^/2m —p 
and = lT?q^/Am — 2/i are the free dispersion relations for fermions and bosons, respec¬ 
tively. Vint(<0) characterizes the binary attractive interaction potential between fermions. 
Additionally, g represents the FR coupling between the closed and the open channel states, 
and uq is the detuning of the boson resonance state from the collision continuum. (We stress 
that the considered approach is actually a realization of a general two-channel model to the 
particular context of Fermion-pair production.) 

Initially, the system is at equilibrium at a finite temperature T. In that situation, a 
sinusoidal modulation of the detuning from the FR is applied. Correspondingly, z/q is re¬ 
placed in Eq. (1) by z/(f) = uq + Asmujpt. It is assumed that which corresponds to 
the pairing interaction resulting from nonresonant processes, is not affected by the applied 
driving field. Through the unitary transformation U{t) = 
the Hamiltonian, transformed as H' = U^HU — ihUW, is converted into 
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( 2 ) 


k,cr 


A 


H' - iiN = X] ( - h—An{!jjpt) 1 + Vm aa+k,t®a-k,;®f-k',^«f+k',t + 


+ hvQ)h\hc^ +{g^ ^q«H+k, 
q V q,k 


q,k,k' 


t4-k,^+h-c. 


(In obtaining the above expression we have made use of the partial result WHU = H, 
which derives from the commutation relation of H with the total number of fermions and 
the cancellation of the introduced time dependence in the interaction term.) (In Ref. [14], 
an alternative approach was implemented by transferring the time variation in u to the 
intermode coupling via a different unitary transformation. Note, that, as shown in Ref. 
[14], the time dependence of the coupling term prevents the one-mode reduction, applicable 
to the undriven dynamics for broad resonances.) An approximate description of the dy¬ 
namics resulting from the modulation can be obtained through the following perturbative 
scheme. The complete Hamiltonian is split as H' — fiN ~ Ho + Hper, where the unperturbed 
Hamiltonian has the form given by Eq. (1), i.e., Hq = H — fiN] and, the perturbation reads 


Hpg^r ^A siu.i^UJpt') ®k,(T®k;,fT • 

k,cr 


(We consider that the modulation amplitude is sufficiently small for the perturbative scheme 
to be valid. It is assumed that the system, initially in the BCS side, stays in that regime 
during the whole process. Hence, the BEC side is not reached and neither is attained the 
unitary limit. Later on, we will precisely define the range of applicability of our approach.) 


A. The zero-order Hamiltonian 


To describe the unperturbed system, we follow the standard HFB approach 18|, |20| . 
Accordingly, we introduce first three mean fields: uq = for the spin density, 

Ao = \Vint \ Xlk (®-k,4.ak,t) for the pairing field, and 0m,o = (^q=o) for the boson field. (We 
take q = 0 as we focus on the condensed molecular field.) Through the incorporation of 
those mean fields, the zero-order Hamiltonian, which describes the unmodulated system, is 
rewritten in the form 


Hq — 'y ^ hfcOk,a-®k,( 

k,(T 


5Z(^o4w-k,t + h-C-), 
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which corresponds to an effective BCS model with mode energy I 4 = + Vint^o and gap 

Aq = Aq — g(j)m,o- The mean-field description inclndes also the equation for the evolution of 
the boson mode, namely, 


^ _ 2/i)0^o _^Ao. 


dt 


\V^ 


( 3 ) 


int 


Hq is straightforwardly diagonalized. By applying the Bogoliubov transformation (BT) 
defined by the fermionic operators Ckq = cos^^Ok^ — sin^^alj^j^ and = sin6'fcak,| + 
cos6*fca^,^j^, where 9^ is given by tan(26'fc) = Aq /14 ISl, l20|, we find 


Ho = '^ Efc,o(Ck_^Ck,t + Ck,;Ck,i) + constant. 

k 

The operator (ck,t) creates (annihilates) a quasi-particle excitation with momentum k 
and spin The associated excitation energies are 


Ek,o = \/V^ + Al= J {tt^k'^/2m - n + VintnoY + A^ 


( 4 ) 


The BCS state |Tbc' 5 ) is the effective vacuum state of this Hamiltonian, i.e., 
Ck,t|TBC 5 ) = CkqlTscs) = 0; in the previous representation, it is given by iTscs) = 
rik (^cos6*fc -|-sin ^j |0) 12|. Note that the excitation gap Aq combines the mean 

pairing field Aq and the equilibrium molecular field 4>rn,o- Aq is obtained from the BCS 
equation 

Ao = ^|2/'»(0) - 1| siti(2ft), (6) 

k 

where {/^'^(O)} are the initial populations of the quasiparticle states, which are given by 
the Fermi distribution function, i.e., /fc'^(O) = 1/(1 -|-since thermal equilibrium is 
assumed for the system before the application of the magnetic modulation. Moreover, (j)rn,o 
is obtained as the stationary solution to Eq. (3), namely, 

_ I/Aq 

0m ,0 — I ,r, -V- (6) 

|4/nt| (2/i — ZZq) 

(Note that a selfconsistent procedure is required to obtain Aq and 4>m,o-) 
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B. The perturbation 


We turn now to analyze the effect of the perturbation. The introduction of in 
the HFB approach implies dealing with changes in the mean helds, which now become 
n{t) = no + Sn(t), A(t) = Aq + SA(t), and = 0m,o + S4>m(t). Our objective is solving 
for the perturbation-induced increments of those helds. In particular, we will focus on 
explaining the damping and delay of SA{t) observed in the experiments. Let us see that 
a hrst-order approximation to that behavior can be obtained simply by incorporating Hp^r 
into the HFB approach dehned by the unperturbed mean helds. Specihcally, we apply the 
previously dehned BT to the perturbation Hamiltonian, which, as a result, is written as 
Hper Hcoup-j where 


Hevar = -^ Sin{ujpt) COs(26'fc) (^4 .f-Ck,t + 4 > (7) 

k 

Hcoup = -^sin(a;pt) ^sin(26'fc) + h-c.) . (8) 

k 

From the forms of H^var and Hcoup-, two preliminary general conclusions on the ehect of the 
magnetic modulation can be drawn. First, Hevar leads to a time variation of the quasiparticle 
energies, which become Ei.{t) = q -|- SE^it), where 

hA 

6Ek{t) = —— sin(ojpt) cos(26'fc) (9) 

In the next section, we will see that the disequilibrium induced by this term and the sub¬ 
sequent relaxation of the condensate are at the origin the detected delayed response of the 
system to the driving. Second, Hcoup represents modulation-induced interactions between 
the vacuum state and a doubly-excited state. Importantly, these coupling terms, which 
oscillate with the external frequency cUp, are relevant only when they can induce an effec¬ 
tive resonance between the BCS state and the two-excitation conhguration, i.e., only when 
hup > 2Ao. The resulting heating effects will be analyzed in Section IV. 
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III. THE EFFECT OF THE QUASIPARTICLE-ENERGY VARIATION 


In the regime defined by tlujp < 2Ao, the interaction terms given by Hcoup can be dis¬ 
carded, and, consequently, the perturbation Hamiltonian Hper can be approximated as 
Hevar = Xlk('^k + '^k^.'^kq)- Heuce, the driven Hamiltonian is still diagonal in the 

representation of the quasiparticle states of the unmodulated system. We will see that, 
although Hevar simply leads to the variation of the quasi-particle energies, the consequent 
effect on the gap dynamics can be quite complex; in fact, its analysis will require the gen¬ 
eralization of our model. (Here, we will follow a treatment alternative to that presented in 
Ref. [14], which will allow us to simplify the characterization of the basic physics of the 
delay.) 

The perturbation forces the system out of equilibrium; the initial thermal populations 
associated with the unmodulated energies do not fit the Fermi distribution = 1/(1 + 

^Eu(t)/kBT'^ for the actual (time-varying) energies. Indeed, the gap equation now reads 

A(t) = ^ ^ [2/fc(t) - 1] sin(20fc), (10) 

k 

where {fk{t)} are the (changing) populations. To describe the dynamics, we must deal 
with the relaxation of the populations towards equilibrium, which implies extending the 
current Hamiltonian description. Note that a selfconsistent approach is needed. The evolv¬ 
ing {Ek{t)} affect the relaxation of the populations {fk{t)} by modifying the equilibrium 
distribution {//"^(t)}. In turn, the variation of A(t) changes the global Hamiltonian, and, 
in particular, can alter the quasiparticle energies. As previously stated, in the simplified 
description considered here, we neglect corrections to the quasi-particle energies due to 
changes in the mean fields: the form given by Eq. (9) is assumed to permanently apply. 
(See Ref. [14] for an analysis of higher-order effects.) We will see that this simplification 
retains the system components responsible for the emergence of the features observed in the 
experiments. 

A. The relaxation mechanism 

The mechanism for thermalization can be assumed to be based on collisions between 
excited particles. Here, instead of tackling a detailed analysis of the dependence of the 
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relaxation on the system characteristics, we will focus on general aspects of its role in the 
condensate dynamics. Accordingly, we consider that the evolution of the populations is 
governed by the generic equation 22j 

% = ~ IA(‘) - A'’(i)l. (11) 

where 1/r/ represents the effective thermalization rate. No restrictions on the magnitude of 
Tf are assumed. A similar relaxation mechanism was considered in Ref. [22] in the context of 
nonequilibrium superconductivity. Central to this mechanism is the idea that the relaxation 
is activated by the distance from the actual populations to those corresponding to the 
equilibrium, which are, in turn, changing as the quasiparticles energies are being modified 
by the driving. We have assumed a first (compact) form of characterizing that process 
with the introduction of the effective thermalization rate. Note that, since the system is 
continuously forced out of equilibrium by the driving field, i.e., the are permanently 

changing, the relaxation mechanism is always activated. The nondirect following to the 
modulation observed in the experiments can be anticipated to be rooted in finite values of 

Eq. (11) is an inhomogeneous linear differential equation, which is exactly solved to give 

Mt} = (^MO) - . 

Furthermore, through integration by parts, we find 

rt Jfsq 

Mt) = (MO) - n-{0)) + Mit) - e-M j e‘'l-i^(t')dt'. (12) 

This expression is simplified by taking /^(O) = /fc'^(O), since the system is at equilibrium 
at f = 0. By combining Eqs. (10) and (12), we obtain the following integral-differential 
equation for the order parameter 


A(f) = ^ 


k >- 


2 fl\t) 


dt 


{t')dt' 


sm(2M)- 


(13) 


We deal now with particular regimes where we can go further in the analytical characteri¬ 
zation of the gap evolution, and, consequently, in the identification of the delay time. 
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B. The response to the modulation at small departure from equilibrium 


Eq. (13) simplifies considerably in the regime defined by ~ A <C T ^ T^, = 

1). (Tc is the temperature for the BCS transition.) In this range, the approximations 


~ /fc'^(O) + -^6Ek{t) and ^ ~ —^ can be made (22|. In turn, we can write 

HA 

~ = -g^ cos(a;pt) cos(26*fc). Through the incorporation of these approximations 

into Eq. (12), and, subsequent integration, we obtain for the populations 


fk{t) = /r(0) +^7X7^^^ [sin(a;7) - UpTf {cos{upt) - e . (14) 

sic i- + [l^pTf) 

Then, combining this equation with Eq. (13), we find 


An 


= 1 -\-C \e sin (p + sin(a;pt — (^)] , 


where 


(15) 


^ ^ \Vint\ - sin(26*fc) cos(26*fc), 

Ao STc^l + {upTfy ^ ^ ^ 


and 


(16) 


ip = aictan^UpTf). 

Some implications of these results must be stressed: 

(i) The gap evolution incorporates a transitory decay with characteristic time tj and a 
secular oscillatory behavior with frequency Up. The external field is not instantaneously 
followed; associated with the phase shift (^, there is a delay time given hy td = ^ = 
Tf[l + 0 {{ujpTf)^)], which can be interpreted as the condensate relaxation time. No changes 
in the delay are observed at different cycles of the field modulation in agreement with the 
experimental results. Furthermore, the detected invariance of the delay with the exter¬ 
nal frequency can be understood as associated with the small magnitude of the correction 
O {{upTf)^) for the experimental conditions. The complex character of the driving mecha¬ 
nism is apparent in the obtained expression for the amplitude of the oscillatory term, which 
combines external-field parameters and characteristics of the unperturbed system. (The 
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factor ^1^8111(20^.) 008(20^.) pre8ent in Eq. (16) can be 8tandardly evaluated 2^, l23|.) It i8 
worth enipha8izing that the effect of the perturbation 8cale8 with the factor Aj 

(ii) Out8ide the con8idered regime, the intricate interdependence of the gap and the pop- 
ulation8 can imply a complex nonlinear contribution of the population8 to the gap relaxation 


22|. Hence, one can expect that, in a general regime, the delay time can 8igni£cantly differ 


from T^. 

(iii) The con8i8tency of our approach can be te8ted by analyzing the limit8 of 8mall and 
large relaxation time When r/ i8 much 8maller than any other characteri8tic time in the 


proce88, in particular, than the driving period, we hnd that ^ 0 


24| . Then, there i8 no 


delay between the gap evolution and the external held. The prediction8 of the adiabatic 
approximation are con8i8tently reproduced: for a 8udden relaxation, the population8 follow 
adiabatically, (on the relaxation time 8cale), the equilibrium value8 {/^'^(t)} a 880 ciated with 
the time-dependent energie8. The a 880 ciated gap dynamic8 become8 “trivial”: the evolution 
corre8pond8 to a 8equence of equilibrium 8tate8 where time enter8 a8 a parameter. On the 
other hand, for a very large r/, we obtain O —)■ 0, and, therefore, A(f) = Aq: we trivially 
recover that there i8 no change in the gap for time8 much 8maller than the characteri8tic 
time for the evolution of the population8 rj. 


IV. HEATING EFFECTS 

For hbjy > 2Ao, the term H^oup in fhe perturbation become8 relevant: the magnetic held 
can then induce an effective re8onance between the fundamental 8tate and a doubly-excited 
8tate. (See Ref8. [25] and [26] for related work.) To analyze the re8ulting tran8ition, we 
rewrite the coupling Hamiltonian a8 


^coup fR 8in(cc:pt) , 


where 


W = —^ ^8in(20fc) + h.c.) . 

k 


Note that, becau8e of the dependence of the factor 8in(26'fc) = 
energy, the coupling i8 le88 effective a8 the energy grow8. 


Aq /Efc 0 on the qua8iparticle 
(The oppo8ite occur8 in the 
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term H^var, [see Eq. (7)].) The combination of this characteristic with the form of the 
density of states will be shown to determine prominent featnres of the system response. The 
transfer from the BCS state I'hsc's) to the doubly-excited state |/) = cj^qc[^k, 4 . ICBC's) of 
the continuum of quasi-particle states can be evaluated applying Fermi’s Golden Rule. (The 
reverse process, i.e., the decay of doubly-excited states induced by Hcoup, can be neglected: 
given the range of temperatures considered, the population of excited states is always much 
smaller than that of the fundamental state. Also relevant to the lack of symmetry in the 
reverse transition is the continuum structure of the excited states.) As corresponds to a 
sinusoidal perturbation, we have for the transition rate: 




W 




BCS 


^~ Ebcs — 


nh 


—A^ ^ sm^(20t)S(2Ekfi - 


(17) 


The sum in k is standardly converted into an integral: —)■ ^ 2 ^ / 0. [Here, V 

denotes the quantization volume, which disappears in the hnal expression as a scaling with 
is incorporated into the dehnition of the operators introduced in Eq. (1)]. The integral 
is evaluated by changing to the variable q- Indeed, incorporating the density of states 
obtained from the dispersion relation given by Eq. (4), the transition rate is found to be 
given by 






271 \2h‘^J 

{hupl2y 


u, 


P L 


{fihJpj2^^ Aq -|- p, Vint'IS'O 


1/2 


X 


^0 


1 - 1/2 


Q{hjjjp — 2 A, 


0 ), 


(18) 


where A = and Q{x) is the Heaviside step function: the excitation process is acti¬ 

vated only for frequencies equal or larger than the threshold value = 2 A 0 . Note that 
the divergence at the threshold is a consequence of the singularity of the density of states 
in the BCS model at Ek^o = Aq. It is important to take into account that this failure of 
the perturbative scheme does not invalidate the identihcation of the physical mechanism re¬ 
sponsible for heating: the fast decrease of the transition rate with Up allows the applicability 
of the used approach sufficiently far from the threshold. 


11 












The transfer of population from the fundamental state to the excited states implies the 
decrease of A(f). This is apparent from Eq. (10): although the total population is conserved 
in the transition, A(t) diminishes since the factor sin(26'fc) decreases as Ek^ grows. Since the 
excitation is permanently activated by the driving, the resulting damping process is contin¬ 
uous, in agreement with the character of the detected decay of the condensate fraction. The 
global picture of the gap dynamics that emerges from combining this effect with the delayed 
oscillation analyzed in Sec. Ill corresponds to the behavior detected in the experiments. 
Note that, in the applied approach, the mechanisms for delay and damping can be consid¬ 
ered to work in parallel. In this sense, it is worth pointing out that the term of population 
gain for the excited states that, because of heating, should be added to Eq. (11) is irrelevant 
given the small contribution of those high-level populations to the gap dynamics. 

The applicability of the perturbative approach can be assessed from the analysis of the 
dependence of the system output on the modulation parameters. Both, the damping coeffi¬ 
cient 7 and the amplitude of the oscillatory component C, given by Eq. (16), diminish for 
decreasing modulation amplitudes and growing frequencies. The observation of oscillatory 
behavior along with damping requires working with decay times larger than the driving pe¬ 
riod. Given the dependence of 7 on and the requirement > 2Ao for the emergence 
of damping, that situation can occur for sufficiently small modulation amplitudes and large 
frequencies. Then, it seems possible to reproduce that situation in a range of parameters 
where the applicability of the perturbative scheme can be guaranteed. Technical details 
of the measurement of the condensate fraction, which is the magnitude reproduced by our 
model, can be found in Ref. [15]. 

We have assumed that the trapping conditions implemented in the referred experiments 
do not crucially affect the main characteristics of the observed features. Indeed, as our 
uniform description qualitatively reproduces the experimental results, the robustness of the 
identified physical mechanisms against spatial non-uniformities can be conjectured. For 
a smooth external potential t/(r), which corresponds to the practical conditions, a local- 
density approximation can be applied to generalize the uniform picture. Accordingly, trap¬ 
ping can be incorporated in the previous framework by replacing the chemical potential as 
H ^ nif) = n — U{f) (27|. The consequent use of local fields implies less compact results 
but does not affect the basic physics underlying the observed features. One of the effects 
of nonuniformity on heating seems evident: the local fields Ao(rO, and, in 
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turn, Ao(t^, can be expected to smoothen the sharp threshold for the onset of excitations. 
Although the previously derived expressions can be modified by the averaging over the 
distribution along the trap, the former basic picture still applies. 


V. CONCLUDING REMARKS 


The considered variation of the coupled fermion-boson model has been shown to give use¬ 
ful clues to understanding the dynamics of atomic fermion pairs with modulated interaction 
strength. The perturbative scheme set up from the basic BCS approach has allowed isolat¬ 
ing the role of the different elements of the system. The emergence of specific dynamical 
features has been found to depend on the time scales of the system components: the relative 
magnitudes of the driving period, the inverse gap frequency, and the relaxation time deter¬ 
mine the characteristics of the system response. The interest of further experimental work 
on specific aspects of the modulation scheme is evident. Particularly valuable can be the 
experimental characterization of the system response for driving frequencies slightly larger 
than the threshold, where the perturbative approach fails. Also interesting can be checking 
the inhibition of damping for frequencies smaller than 2Ao. 

It is of interest to mention recent related works on alternative approaches to the study 
of similar systems. In this sense, it is pertinent to point out the advances in the charac¬ 
terization of the system parameters which are reported in Ref. [31] and the applications of 


analogue methodology to p-wave interacting Fermi gases 


32j. Also valuable is to establish 


a parallelism between the considered scenario and similar modulation techniques applied in 
bosonic systems. (Ref. [33] presents interesting experimental findings on the production of 
ultracold molecules via a sinusoidal modulation of the magnetic field. Subsequent theoretical 
analysis was presented in Ref. [34].) 

Finally, it is worth pointing out that the applicability of the study is not restricted 
to the field of ultracold atomic gases. The central issue in the analysis, namely, the gap 
dynamics 28|, in particular, the effect of changes in the populations on the evolution of 
the condensate fraction, is relevant to topics ranging from nonequilibrium superconductivity 


22 


. l29| to quenched dynamics in superfluid ^He |30|. In this sense, the analytical character 


of the study can be particularly useful given the difficulty of dealing with out-of-equilibrium 
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situations in those contexts. 


[1] For recent reviews, see S. Giorgini, L.P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 
(2008); I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008); M. Lewenstein, 
A. Sanpera, V. Ahufinger, B. Damski, A. Sen De, and U. Sen, Adv. Phys. 56 , 243 (2007); and 
references therein. 

[2] M. Greiner et ai, Nature (London) 426 , 537 (2003). 

[3] S. Jochim et ai, Science 302 , 2101 (2003). 

[4] M. W. Zwierlein et ai, Phys. Rev. Lett. 91 , 250401 (2003). 

[5] T. Bourdel et al., Phys. Rev. Lett. 93, 050401 (2004). 

[6] M. Bartenstein et al., Phys. Rev. Lett. 92 , 120401 (2004). 

[7] C. A. Regal et al., Phys. Rev. Lett. 92 , 040403 (2004). 

[8] M. W. Zwierlein et al., Phys. Rev. Lett. 92 , 120403 (2004). 

[9] K. E. Strecker et al., Phys. Rev. Lett. 91, 080406 (2003). 

[10] M. Bartenstein et al., Phys. Rev. Lett. 92 , 203201 (2004). 

[11] J. Kinast et ai, Phys. Rev. Lett. 92 , 150402 (2004). 

[12] W. Ketterle and M. W. Zwierlein, Proceedings of the International School of Physics ’’Enrico 
Eermi”, Gourse CLXIV, (lOS Press, Amsterdam, 2008). 

[13] E. Taylor et al., Phys. Rev. A 74 , 063626 (2006); N. Eukushima et al., Phys. Rev. A 75 , 033609 
(2007). 

[14] J. Plata, Europhys. Lett. 87, 50001 (2009). 

[15] M. W. Zwierlein et ai, Phys. Rev. Lett. 94, 180401 (2005). 

[16] M. Holland et al., Phys. Rev. Lett. 87, 120406 (2001). 

[17] E. Timmermans, K. Furuya, P. W. Milonni, and A. K. Kerman, Phys. Lett. A 285, 228 (2001). 

[18] M. L. Chiofalo et al., Phys. Rev. Lett. 88, 090402 (2002). 

[19] Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89 , 130402 (2002). 

[20] Y. Ohashi and A. Griffin, Phys. Rev. A 72 , 013601 (2005). 

[21] Y. Ohashi and A. Griffin, Phys. Rev. A 72 , 063606 (2005). 

[22] M. Tinkham, Introduction to Superconductivity (Dover Publications, New York, 1996). 

[23] M. Marini et al., Eur. Phys. J. B 1 , 151 (1998). 


14 



[24] W. Yi and L. M. Duan, Rev. A 73, 013609 (2006) 

[25] M. Greiner et al., Phys. Rev. Lett. 94, 070403 (2005). 

[26] J. Plata, Phys. Rev. A 74, 013603 (2006). 

[27] L. Salasnich, N. Manini, and A. Parola, Phys. Rev. A 72, 023621 (2005). 

[28] R. G. Scott, F. Dalfovo, L. P. Pitaevskii, and S. Stringari, Phys. Rev. A 86, 053604 (2012). 

[29] J.A. Schmid and G. Schon, J. Low Temp. Phys. 20, 207 (1975). 

[30] G. L. Warner and A. J. Leggett, Phys. Rev B 71, 134514 (2005). 

[31] M. H. Szymahska et al, Phys. Rev. A 72, 013610 (2005). 

[32] L. Austen et al., Phys. Rev. A 87, 023610 (2013). 

[33] S. T. Thomson et al., Phys. Rev. Lett. 95, 190404 (2005). 

[34] T. M. Hanna et al., Phys. Rev. A 75, 013606 (2007). 


15 



